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Abstract 

Secondary radiation of hadrons from a tau pair produced in electron positron collisions may con- 
, stitute an important obstacle for precision measurements of the production cross section and of 

branching ratios. The rate for real and virtual radiation is calculated and various distributions are 
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presented. For Z decays a comprehensive analysis is performed which incorporates real and virtual 
radiation of leptons. The corresponding results are also given for primary electron and muon pairs. 
Compact analytical formulae are presented for entirely leptonic configurations. Measurements of 
Z partial decay rates which eliminate all hadron and lepton radiation are about 0.3% to 0.4% 
lower than totally inclusive measurements, a consequence of the 0(a 2 ) negative virtual corrections 
I/*-) , which are enhanced by the third power of a large logarithm. 

(N 
O 
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1. Introduction 

i \ Since the turn on of LEP several years ago experiments have collected an impressive amount of data, 
' with a combined sample of more than 10 7 hadronic and 10 6 leptonic events. This enormous event 
rate has allowed us to measure the ratio of hadronic versus leptonic events or the leptonic decay rate 
with an accuracy of better than two per mill. Clearly, at this level of precision subtle aspects of 
^ ■ radiative corrections have to be taken into account. One specific example is the Z decay rate into 
a lepton pair: Details of the treatment of additional radiation of a pair of soft leptons or of a soft 
hadronic system have a strong impact on the predicted rate. Formally the rate of these events is 
of order (a/ir) 2 ~ 5 • 1(P 6 . In reality, however, these events are enhanced by the third power of 
a large logarithm, which essentially compensates one factor of a/ir. This positive contribution is 
to a large extent cancelled by the corresponding virtual correction, leaving a small positive shift in 
the inclusive rate. The dominant part of this shift can be incorporated into the corrections of order 
a /it by employing the running QED coupling a(Q 2 ), the remainder is truly negligible. However, in 
an experimental analysis which eliminates events with secondary radiation the large negative virtual 
corrections must be taken into account. 
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With this motivation in mind we present in this paper a comprehensive study of virtual and 
real leptonic and hadronic radiation. Particular emphasis is put on the radiation of hadrons from a 
primary r lepton pair. This is the most difficult reaction to evaluate, since three independent scales 
s, m T and rrihad enter the calculation. The results for the radiation of hadrons from an electron or 
muon pair are significantly simpler: The lepton mass is far smaller than the threshold for hadron 
production and can be safely set to zero throughout the calculation. Results for this latter case have 
been given in [jij for virtual corrections in the high energy approximation. The corresponding case 
for real final state radiation with a massless primary lepton has been solved in analytic form in ||, [|] . 
The virtual and real radiation of a lepton pair off a massive primary lepton has been evaluated for 
arbitrary s in [Q] for the special case that the mass of the secondary lepton is far smaller than all 
other scales of the problem. In this paper we complete this topic by considering radiation of hadrons 
from a primary r pair. Since large data samples of r pairs are also collected at lower energies, e.g. 
at around 10 GeV, this case also is scrutinized. However, it is observed that radiation is far less 
important in this low energy region, an obvious consequence of the less prominent role of the "large" 
logarithms. Hadronic radiation off primary r pairs also plays a special role from a purely practical 
point of view: t leptons decay into hadrons (plus a neutrino) and hence events with additional hadron 
radiation might well be considered experimentally as hadronic final states, since they have larger 
hadronic multiplicity and the invariant mass of the hadronic subsystems exceeds m T . By the same 
token, the measurement of r decays into large multiplicity final states could in principle be influenced 
by inappropriate assignment of the radiative events. Last but not least, signals of "new" physics could 
be faked by a misinterpretation of these Standard Model reactions. 

The final state radiation discussed in this paper is complementary to the initial state radiation 
treated in |j| . No interference terms are present between radiation from the initial and final state as 
long as the (dominant) QED reactions are considered and axial vector couplings from neutral current 
reactions are ignored. Final state radiation of lepton (and quark) pairs has also been treated in || in 
the context of a Monte Carlo Program. For hadron radiation this approach turns out to be inadequate, 
since the reaction is dominated by low mass hadronic systems. For lepton radiation our results are 
lower than those of || by about 15%. An important ingredient for the evaluation of real and virtual 
hadronic radiation is a proper description of R(s), which has to be taken from experiment. In our 
approach we employ the most recent parametrization of data taken from |7| . 

The plan of this paper is as follows: In section 2 the "master formula" for calculation of the 
total rate is presented, which is at the same time a convenient starting point for a variety of differential 
distributions. The case of real radiation with primary r- leptons and secondary hadrons is discussed 
in some detail, considering high energies (y^ = Mz) as well as yfs = 10 GeV (characteristic for a B 
factory) and 4.2 GeV (relevant for a r-charm factory). Predictions for the total rate are presented, 
as well as distributions which are differential with respect to the mass of the hadronic system and its 
energy. The corresponding virtual radiation is calculated in section 3 for the same cm. energies and 
at the r pair threshold. Corrections are given for the Dirac and Pauli form factor, as well as for the 
total rate. Section 4 is concerned with a detailed discussion of Z decays. Results for electron, muon 
and tau as primary leptons are presented both for hadronic and leptonic radiation. The corresponding 
predictions for virtual radiation are also collected. It is demonstrated that these virtual corrections 
reduce the exclusive leptonic rate of the Z boson by several per mill. Section 5 contains a summary 
and our conclusions. 
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Figure 1: Typical diagrams for real and virtual emission of hadrons in r pair production. 



2. Real radiation 



The calculation of the ratio o- T + T - had /a p t (where a p t = 47ra 2 /3s) can be reduced to the numerical 
evaluation of the two dimensional integral 

2 (v/i-2m T ) 2 / 

R T+T - had = °-2llZ^ = \(°) I ^ R had {s') F(s'/s) , (1) 
o-pt 3 \ttJ J s' 
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F{s' / s) can be interpreted as the normalized cross section for r pair production with the additional 
emission of a vector boson of mass \fs' '. The integral in (|l[) sums up all hadronic contributions of mass 
yfs' weighted by Rhad(s')- The numerical treatment of this integral allows us to use experimental input 
for the function Rhad(s') and avoids the poor approximations of a parton model inspired calculation. 
In contrast to [II — there we solved the integrals analytically in the case where the "massive" photon 
splits into a pair of fermions with mass far smaller than m T — F{s' / s) is calculated numerically. We 
would like to emphasize that (|l]) contains no restriction on the mass values and is thus also applicable 
to completely different kinematical situations. In the limit of vanishing external mass m T = we 
recover eq. (2) of |J for F(s'/s). 

The physical situation described by ([!]) consists typically of a r pair with large invariant mass 
and a hadronic state with low invariant mass (see also the discussion of our results below). The 
hadrons are emitted relatively collinear with one of the r leptons. The reverse configuration, where 
a r pair is radiated off the corresponding quark via a virtual photon could be treated with a similar 
technique. This final state would consist mainly of quite soft and collinear virtual photons leading to 
a r pair of low invariant mass and is therefore easily distinguished from the situation we focus on in 
this letter. As far as the total rate is concerned, the interference term between these two amplitudes 
vanishes if the production is induced by the (electromagnetic or neutral) vector current. For axial 
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vector induced contributions the interference term is different from zero. However, this term does not 
exhibit a mass singularity in the limit of vanishing quark or leptonic mass and thus is not enhanced 
by a large logarithm. It can be neglected in all cases of interest. In the limit of vanishing fermion 
mass this interference term has been calculated in 

For most of this paper we will discuss the total rate for secondary radiation. However, eq.((l]) 
contains implicitly also information on single and double differential distributions. The integration 
variable s' denotes the squared hadronic mass. The square of the invariant mass of the primary t + t~ 
system, on the other hand, is given by y ■ s, which can be easily related to the energy Ehad of the 
hadronic system. The single differential distribution dR T + T - had /ds' is trivially obtained from eq. ([]]). 
The energy distribution of the radiated hadronic system is given by 
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with 



sy = s + s' - 2E hadx /s and s' min = max (Am 2 . - s + 2E had y/s , 4m£J . (3) 

These formulae will be useful below. 

In the numerical evaluation of (|l]) we use experimental data for Rhad{s') as presented in ]7j[]. A 
convenient approach adapted to the experimental procedure is to decompose the hadronic final states 
into continuum and resonance contributions. Following |7||, the two pion contribution (=p) is taken 
into account together with the continuum by a direct integration of data points. The most important 
narrow resonances are added separately using the narrow width approximation^, which is accurate 
enough for our purposes. More details about Rhad(s') are given in the appendix. 

In Table |l] we display the continuum plus p and the resonance contributions separately for 
different center of mass energies y/s < Mz- The contribution from a virtual Z as intermediate state 
remains unimportant up to y/s ~ Mz- They are strongly suppressed for obvious kinematical reasons 
as can be easily seen in Table |] (see also discussion below). This in turn justifies neglecting radiation 
through the virtual Z. In Fig. ^ R T +r-had ls plotted for the range 2m T < \fs < Mz- The rise of 
-had an d R t +t-i+i~ f° r increasing y/s is easily understood from the leading term in the high 
energy limit proportional to the third power of the logarithm of s, i.e. ln 2 (m 2 iad /s)ln(m 2 /s), rrihad 
being a characteristic hadronic mass scale of 0(300 MeV). (For leptonic radiation these logarithms 
can be calculated explicitly, see eq. @.) 

As an interesting exercise one may apply ((!]) to T + r~qq production where the quarks are again 
radiated off the tau leptons leaving r pairs of large invariant massf]. This leads to a parton model 
like prediction of R T + T -had- Adopting the quark mass values of Q, viz. m u = m d = m s = 300 MeV, 
m c = 1.5 GeV and m& = 4.5 GeV and choosing y/s = 10 GeV / Mz one would predict 

1 We thank F. Jegerlehner and S. Eidelman for providing their compilation of data and discussing its proper treatment. 

2 Masses and electronic widths are taken from l|9|. 



3 In this case one just has to use R(s') = N c Q 2 . (1 + —p 3 -) y 1 — -^t 2 - in the numerical integration. 
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Table 1: Contributions to R T + T ~had an d values of R T + T -i+i- (I = e, (j, ,t) for center of mass energies 
\fs = 4.2 GeV/10 GeV /Mz- y 's' max is the upper limit of the phase space integral ^). 




10 20 30 40 50 60 70 80 90 

Vs [GeV] 



Figure 2: Normalized cross section R T + T -had- 
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Figure 3: Dependence of the normalized cross section R T + T -had from a cut s' cut for yfs = 10 GeV and 
Vs = M z . 



Rr+r-uu = 4R T +T~dd =4Rt+t-ss = 1-1 -lO" 5 / 2.1 -10- 4 , 

R T+T - c5 = 1.2-10" 7 /5.1-10~ 5 , 
# T+T - fe 5 = /2.6-10- 6 , 

^r+r-had = 1.7 • 10" B / 3.7 • 10~ 4 . (4) 

It is obvious that this model underestimates the production rate by about a factor of two. (The same 
phenomenon has also been observed for hadron radiation off muons and electrons H|.) For energies 
near the t + t~ production threshold a parton model like prediction is completely impossible without 
artificial tuning of the masses of the light quark flavours. 

For completeness we also list in Table [l] the rate for the radiation of lepton (e, fi, r) pairs. The 
prediction for the radiation of e + e~ and is based on the analytical formula given in Q, the one 

for t + t~t + t~ is obtained through numerical integration. An adequate approximation (better than 
1%) for R T+T - T+T - at yfs = Mz is presented in eq. (HJ) at the end of section 4. 

Table [l] and Fig. Q demonstrate that hadronic radiation may become important for precision 
measurements of r branching ratios on top of the Z. Radiation of e + e~ or may even be of 

relevance in the 10 GeV region, e.g. at a B meson factory. 

It is evident that experimental sensitivity and the probability to assign r + r~ events with 
hadronic radiation to exclusive or inclusive r pair production depends strongly on the invariant mass 
of the hadronic system. As mentioned above the integral (jl|) is in fact dominated by small s' . In the nu- 
merical approach this is easily shown by replacing the upper limit of the integral in (|I|) by a value s' cut < 
s'max = {yfs-2 m T ) 2 . In Fig. | we display our results for R T+T - had (s, s' cut ) = a T+T - had (s, s' cut )/a pt for 
yfs = 10 GeV/M z . Even for Z decays 86% of the integral originates from \fs' < 5 GeV and 58% from 
\fp < m T . In Fig. || we present the differential distributions dR T + T ~ had /ds' and dR r + T - had /dEf iad for 
yfs = 10 GeV / Mz using eqs.(^j) and (^). Contributions from the narrow resonances are not taken into 
account for the distributions dR T + T ~h ad /ds' but can easily be read off from Table [l|. As expected the 
distributions are peaked at small s' (already evident from Fig. ^) and at fairly small Eh ad illustrating 
that hadronic radiation is clearly dominated by the p-resonance. It is not inconceivable that these 
events might affect precision measurements of multipion decay modes of the r. In view of possible 
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Figure 5: Differential distributions for real emission of e + e in r pair production as described in the 
text. 

experimental applications we also present in Fig. || the differential distributions dR T + T - e + e - /ds' and 
dR T + T - e + e - /dE e + e - for the radiation of electron-positron pairs at yfs = 10 GeV and at Mz- These 
are evidently also peaked at extremely small values of s' and E e + e - . 

If one would simply eliminate all events with secondary radiation of hadrons and leptons, one 
would arrive at a measurement of T{Z — > t + t~) significantly below (about 0.4%) the generally quoted 
inclusive value. Both real and virtual radiation are enhanced by the square of the large logarithm 
ln(M|/m? 0< /) times ln(M|/m?) . Because the coefficients of these logarithms only differ by sign, real 
and virtual corrections cancel to a large extent. This leads to a small contribution to the inclusive 
decay rate, as can be seen in the next section. Measurements of exclusive r pair production, however, 
are strongly affected. With this application in mind we proceed to the calculation of the virtual 
corrections. 
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Table 2: Virtual corrections due to hadronic and leptonic vacuum polarisation for center of mass 
energies yfs = 2m r /4.2 GeV/10 GeV/Mz for form factors and the resulting cross section as described 
in the text. 



3. Virtual radiation of hadrons 

Virtual corrections to r pair production affect both Dirac and Pauli form factors, F\ and F%. Leptonic 
corrections have been calculated in analytic form in j|] for the cases that the masses of the radiated 
leptons are much smaller than or equal to the r mass. The effect of virtual radiation of hadrons 
via hadronic vacuum polarisation of one exchanged photon, as depicted in Fig.[l]b, is given by an 
expression similar to eq. (|l|). For virtual radiation the phase space integration has to be extended 
to infinity corresponding to the dispersion integration of the absorptive part of the hadronic vacuum 
polarisation. The function F(s',s), now representing the normalized cross section for the production 
of a r pair with an exchanged vector boson of mass x/s 7 , may be easily calculated analytically and 
reads 



F(s',s) 

where 



a^ 2 r 



7T 



lu(3-lu 2 ) Reims' ,s) +ReF 2 {s',s) + to 6 ReF 2 {s', s) 



(5) 



uj = Jl-Am 2 T /s. (6) 



\ 2 f ds' 

J J —R ha d(s')F( S ',s) (7) 



The functions Reiq^ are given in ||] (eqs. (15) and (16)) and denote the real parts of the electromag- 
netic form factors arising from the exchange of a vector boson of mass \fs'. The integral 

R « rt _ i (<*V r ds' 

is calculated numerically. A similar formula applies to the form factors individually. 

Virtual corrections from different hadronic intermediate states cannot be separated experimen- 
tally and contribute in a coherent manner. Their combined effect on the form factors and the rate 
are displayed in Table ^ for four characteristic energies. For completeness we also display the con- 
tributions from leptonic vacuum polarisation. Contributions from virtual electron pairs are obtained 
with analytical formulae derived in the limit m e <C m T (see In the same reference an analytic 
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Figure 6: Contribution to e + e~ — * t + t~ due to virtual radiation of hadrons. 

formula for virtual r pairs is also at hand. It is evident from Table ^ that these corrections must 
be taken into consideration in precision measurements of partial rates. Close to threshold the virtual 
corrections give a positive contribution to R as is illustrated in Fig. ||, where R^+ t T - had is plotted 
for the range 2m T < \/s < Mz- This is a consequence of the Coulomb attraction mediated by the 
exchange of a vector boson. Far above threshold the corrections are negative. They lead to a reduction 
of the exclusive decay rate of Z — > t + t~ by nearly three per mill, exceeding the 0(a) correction of 
3/4(a/7r) ~ 1.7 • 1CP 3 . As mentioned above these negative corrections cancel for the most part the 
positive corrections from real emission, leaving a small positive overall correction. This remainder 
can to a large extent be absorbed into a running electromagnetic coupling constant a(s) as discussed 
in@]. 

4. Corrections to Z decays into and e + e~ 

In jl|, ^] analytical and numerical results were presented for real and virtual radiation of massive leptons 
or hadrons from the primary production of massless leptons. For completeness it seems appropriate to 
present in this paper a brief update of the hadronic radiation in Z decays to /^ + /U~ and e + e~ together 
with a complete table for the radiation of leptons. Real and virtual radiation are presented in Tables [| 
and |H respectively. It is evident that exclusive determinations of Z — > e + e~ or will become 

sensitive towards these corrections. 

The results for R e + e -had an d R^+^-had are ~6% larger than our old value ||. Half of this 
effect is due to the updated description of the continuum and 3% originate from the inclusion of 
the effective electromagnetic coupling constant a (see also the appendix). Comparing our result for 
Rg+e-had (taking y/s = 93 GeV) with the corresponding one in [jij the difference amounts to only 
1%. Our results for the radiation of leptons are about 15% lower than those of ||, essentially a 
consequence of the different choice for a. There the running electric coupling constant a(M^) was 
chosen. However, radiation of leptons, in particular of electrons, is dominated by small s'. Therefore 
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Table 4: Virtual radiation R^Lary, secondary f or s = M 2 z . 



the low energy scale adopted in this work is more appropriate. 

Below we also list the high energy expansions (neglecting terms of order my 2 /s) of Rf 1 f 1 f 2 f 2 
for the cases m\ S> ni2, mi = ni2 and mi <C TO2, where mi is the mass of the primary 



produced fermion fi and m<i the mass of the secondary emitted fermion fi- These formulae provide 
excellent approximations for \fs~ = Mz (better than one per mill in all cases except for fi = r, where 
the relative difference between complete and approximate result amounts to nearly 1%) and therefore 
can directly be applied to the corresponding leptonic Z partial widths. 



ft m,i>m 2 
/1/1/2/2 



ft m_l=m_2 
/1/1/2/2 



ft m\<^m 2 
/1/1/2/2 



-In^ifln^i 

s 6 V s 



+ 



+ 1 
13 



In ™1 (- In 

n s V6 n s 18 s ' 3 sv " ; 36 



2 rnf_13 mf 4 53; 
s 18 s + 3 U ' Sfi 

2 . , \ , m\ 

s 



: 36 s \ 1 08 3 s v ; / .<? 



18 s 36 s 
5 . 32 w , 833 3 mo 9 

+ i«») + T« 2 »-5iH-i^^ 

73 



/«\ 2 r L i n 3 ZZll _ 1^ i n 2 2/73 

W I 18 n s 36 n s + 3 V 18 + ^ 



18 s 

, 11 w . 1829 
+ C(3) + T C(2)- W 



1 . o mn 19 , 9 

l — — In 2 — + 



In 

18 



36 



+ |C(3) + f C(2) 



s 

2123 
324 



-| + c, 2 , 



In^t 



In^t 



(8) 



(9) 



(10) 



11 



R 



hhhh ^ 



In 2 — - In — + 1 

s 6 V s 



, m?/ 1 , o m? 13 , m? 4 , 11\ 
+ ln^ - ln 2 ^ + - ln^--C 2 + — 
s V 6 s 18 s 3 9/ 



1 , 3 m 13. 2 m? /133 2 ..A , m\ 

H In 01 — In 2 — + h - C(2) In — 

18 s 36 s V108 3 v 7 s 



§«*> 



32 w , 67 3 mo 9 

— C 2) + — + -vr 2 

9 sv ; 27 4 mi 



(11) 



mi=rri2,virt 

hhhh 



R 



rai<^rn,2,virt 

hhhh 



\ 1 , o ml 19 , m\ 2 /265 A , m| 11 A/ , 383 1, , 
-ln 3 ^ + -ln 2 ^ + - (2 In-^--C(2) + — - 12) 

.18 s 36 s 3 V 72 / s 3 54 



18 ~ s ' 36 

1 , ,m| 19 , 2 m| 2 /265 > . A , m\ 
§C(3)-f C(2) + 



s 3 V 72 
3355 



648 



(13) 



Results for virtual corrections were first obtained in ||], [10] and [gj] for the cases mi 3> m2,mi = m2 
and mi <C m2, respectively. Those for real radiation are taken from for mi 3> m,2 and from B ^] 
for mi <C m2- The formula for real radiation with mi = m2 is new. The leading cubic logarithms for 
all these cases can also be found in [11]. Note that in (JsJ) and (|il] ) linear mass corrections ~ m2/mi 



are included. They significantly improve the quality of the approximations, especially in the case 
mi = m T , m2 = m«. Neglecting these linear mass corrections the value of say RT^^ Tt for the 
energy y/s = 20 GeV differs from the value based on numerical integration by 4%, while including 
them leads to a difference of only three per mill. Of course, for higher energies the latter difference 
decreases and amounts to eight per mill vs. half a per mill for y/s = Mz- There are no linear mass 
corrections to the other high energy approximations. It is an interesting fact that the complete results 
for these linear mass corrections, valid for all energies above threshold, are remarkably simple: 
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(14) 
(15) 



with uj defined as in (||). They are in very close relation to the non-analytic linear vector boson mass 
terms of the corresponding one loop results denoted by F, eqs. (||) and (||). The difference between 
fll4| ) and ( |l5| ) can be reproduced by the analogous linear mass correction in the two loop relation 
between pole and running masses, see [12]. 

Again, concerning the total Z decay rate, it is obvious from eq. (|8|)a'0) that by adding the 
corresponding real and virtual contributions the dominant logarithmic terms cancel, leaving the result 



r> _ i nvirt 

u hhhh + K hhhh 



iln4 + C(3)-^ 
4 m-2 8 



(16) 



which is universal for all three mass constellations in the high energy limit. As mentioned in the 
previous section the remaining logarithm can be absorbed into the running coupling constant a(s), 
leaving a tiny term ~ (a /vr) 2 without logarithm. 

For completeness we also present the prediction for the radiation of hadrons from primary 
electrons and muons. In this case the result can be expressed in terms of the moments 



pi dx ln n x / 

Rn= j- [Rhad(4m 2 /x) - Rhad(oO 

Jo x n\ \ 

leading to the following formulae for real 0] and for virtual Q radiation (I = e,fx) 
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(19) 



The evaluation of (|i~7|) gives 
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Rq 

Ri 
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4.31. 
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-48.00 



(20) 



where we used the dressed Rhad{s') including the narrow resonances and vs' = 100 GeV as the 
integration cutoff. The values of the moments individually depend on the choice of this cutoff (and 
correspondingly Rhad(oo)), but we have checked that the results of fll8| ) and (|i~9|) are identical for all 
cutoff energies above 40 GeV. For the derivation of eqs. ([lS|) and (|l9|) the approximation mi <C m n 
was essential. The relative difference between the asymptotic formulae and the complete numerical 
results is at most 0.2% for ^/s = Mz- For lower energies, of course, the difference increases. Similar 
formulae for radiation off the primary tau pair provide a poor approximation to the true answer for 
real and virtual radiation individually, but are adequate for the sum, which is identical to the electron 
and muon case (see also Tables || and ||). For the sum of real and virtual radiation of leptons and 
hadrons at y/s = Mz, which is valid for all three lepton channels f\ = e,fj,,r, one arrives at 
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where the usual definition of the running electromagnetic coupling constant 
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a(s) 



- = n( a ) 
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(21) 



(22) 



was employed. Here the hadronic contributions to the photon vacuum polarisation have also been 
parametrized by the moments defined in fll7|): 



n( s ) ^ - J- 

O 7T 
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Ami 
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(23) 



The differential distribution dRi+i- had /ds' had for the case of massless leptons is given in closed 
form (see eqs. (1) and (2) of M), dRi+i-^/dE^, of course, has to be calculated numerically. The 
predicted distributions for the radiation of hadrons and of e + e~ pairs off primary produced muon 
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Figure 7: Differential distributions for real emission of hadrons and of e + e pairs in fj, pair production 
for v/i = M z . 



14 



pairs are shown[] in Fig. 0. Again the dominant configuration consists of an with invariant mass 

fairly close to Mz and a soft hadronic or e + e~ system. 

5. Summary 

A comprehensive treatment of the radiation of hadronic final states from a primary r pair is presented, 
which is applicable over the full energy range of interest in the foreseeable future. In the low energy 
region real as well as virtual radiation is small, typically at the level of 10 -4 compared to the primary 
reaction. The situation is entirely different for a cm. energy around 90 GeV. In this case important 
differences arise between "inclusive" and "exclusive" measurements of the Z — > t + t~ rate, in particular 
if radiation of the various hadronic and fermionic channels is collected, or eliminated by experimental 
cuts. We have presented differential distributions with respect to the invariant mass and the energy 
of the hadronic system, and for comparison the same quantities for the radiation of electron-positron 
pairs. It is evident that all distributions are peaked at small masses and energies respectively. The 
calculation of the rate for hadronic radiation must necessarily rely on a numerical integration of 
experimental data. Leptonic radiation, on the other hand, can be calculated analytically for the 
mass assignments of interest. Handy formulae are presented for different mass assignments which are 
applicable in the high energy limit, and which are valid up to terms of order mf ept /s. They may be 
useful also for other applications. 
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Appendix: R had {s') 

In our calculation we used experimental data for Rhad(s') as presented in a recent publication by 
F. Jegerlehner and S. Eidelman For a detailed discussion we refer the reader to this work. 

In the region close to the 7r + 7r~ threshold where no experimental data are available we use the 
chiral expansion of the pion form factor (eq. (18) of [0]). The continuum- and p-contributions are 
parametrized by weighted averages of experimental data from various experiments. From 318 MeV up 
to 810 MeV these are given in terms of the pion form factor l-F^I 2 - The continuum from 810 MeV up 
to 40 GeV (above which we use the prediction of perturbation theory) is well described by weighted 
averages of Rhad( s )- At this point two comments are in order: The higher resonances -0(4040), "0(4160) 
and -0(4415) are included in the weighted averages of the continuum data (as can be seen in Fig. 7 
of JjJ) whereas the resonances u, J/ip(lS), -0(25), -0(3770) and the T-family (six resonances) are 
treated in the narrow width approximation as described above. 

All experimental values of Rhad (and l-F^I 2 ) are given as "undressed" (= lowest order in QED) 
quantities. In order to resum all additional self-energy insertions in the photon propagator we use for 

4 The figure for dR l J rl - had / ds' had does not incfude the narrow resonances. The missing contributions to Ri+i-u ad are 
id: 2.33 • 10" 5 , 0: 3.56 ■ 10" 5 , J/ip: 2.95 ■ 10" 5 and T: 8.33 • 10" 7 . 
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the integration of (|T]), (Q) and (^) the "dressed" quantity 

RtT ed (s) = (^) 2 R had (s) , a(s) = . (24) 

V a ) 1 — Aa 

In the low energy region up to v? = 3 GeV most experiments only used the leptonic contribution of 
Aa to extract one-particle-irreducible quantities from the measured cross-sections. For that reason 
we take Aa = Acti ep in this region. For all other energies 

Aa = Aa iep + Aa had . 

As mentioned above these corrections led to an increase of the result by a factor of ~ 1.03. We 
also calculated the effect of Aa; ep (s) on our results, comparing the full one loop expression of the 
fermionic part of the photon self energy with the leading logarithmic behaviour (as it appears in the 
renormalization group "running a"). The difference of ~0.5% may be considered as a crude estimate 
of other effects not taken into account via the implicit resummation by using a(s) and is smaller than 
the claimed accuracy of our calculation. 
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